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Abstract. We derive asymptotics for the moments as well as the weak limit of the height distribution of water- 
melons with p branches with wall. This generalises a famous result of de Bruijn, Knuth and Rice [4] on the average 
height of planted plane trees, and results by Fulmek [9] and Katori et al. [15] on the expected value, respectively 
higher moments, of the height distribution of watermelons with two branches. 

The asymptotics for the moments depend on the analytic behaviour of certain multidimensional Dirichlet series. 
In order to obtain this information we prove a reciprocity relation satisfied by the derivatives of one of Jacobi's 
^ , theta functions, which generalises the well known reciprocity law for Jacobi's theta functions. 
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1. Introduction 



The model of vicious walkers was introduced by Fisher [7]. He gave a number of applications in physics, such 
as modelling wetting and melting processes. In general, the model of vicious walkers is concerned with p random 
walkers on a d-dimensional lattice. In the lock step model, at each time step all of the walkers move one step in 
any of the allowed directions, such that at no time any two random walkers share the same lattice point. 

• A configuration that attracted much interest amongst mathematical physicists and combinatorialists is the 
watermelon configuration, which is a special case of the two dimensional vicious walker model. See Figure [T] for an 
example of a watermelon, where, for the moment, the broken line labelled 13 should be ignored. This configuration 
can be studied with or without presence of an impenetrable wall, and with or without deviation. We proceed 
with a description of p '-watermelons of length 2n with wall (without deviation), which is the model underlying this 
paper. Consider the lattice in R 2 spanned by the two vectors (1,1) and (1,-1). At time zero the walkers are 
located at the points (0, 0), (0, 2), ... , (0, 2p — 2). The allowed directions for the walkers are given by the vectors 

j (1, 1) and (1, —1). Further, the horizontal line y = is an impenetrable wall, that is, no walker is allowed to cross 

• the x-axis. The walkers may now simultaneously move one step in one of the allowed directions, but such that at 
no time two walkers share the same place. Additionally we demand that after 2n steps all walkers are located at 

(M* (2n,0),(2n,2),...,(2n,2p-2). 

Tracing the paths of the vicious walkers through the lattice we obtain a set of non-intersecting lattice paths with 
, steps in the set {(1,1), (1,-1)}. In the case of watermelons without deviation, the i-th lattice path, also called 
T". ' i-th branch, starts at (0, 2%) and ends at (2n, 2i). Further, it is seen that the bottom most path is a Dyck path, so 
^ ■ that the 1-watermelons with wall correspond to Dyck paths. 

k> ' Since Fisher's introduction [7] of the vicious walkers model numerous papers on this subject have appeared. 

While early results mostly analyse vicious walkers in a continuum limit, there are nowadays many results for 
C$ ' certain configurations directly based on the lattice path description above. For example, Guttmann, Owczarek 
and Viennot [13] related the star and watermelon configurations to the theory of Young tableaux and integer 
partitions, and re-derived results for the total number of stars and watermelons without wall. Later, Krattenthaler, 
Guttmann and Viennot [18] proved new, exact as well as asymptotic, results for the number of certain vicious 
walkers with wall. Recently, Krattenthaler [17] analysed the number of contacts of the bottom most walker in the 
case of watermelons with wall, continuing earlier work by Brak, Essam and Owczarek [22] . 

In 2003, Bonichon and Mosbah [2] presented an algorithm for uniform random generation of watermelons, which 
is based on the counting results by Krattenthaler, Guttmann and Viennot [18] (see Theorems 1 and 6 therein). 
Amongst other things they used their generator for obtaining experimental results on the height of watermelons. 
Here, the height of a watermelon is defined as the smallest number h such that the upper most branch does not 
cross the horizontal line y = h. See Figure [T] for an example with four branches and height 13. 
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FIGURE 1. Example of a 4- watermelon of length 18 with wall and height 13 

As already mentioned, watermelons with one branch are simply Dyck paths. It is well-known that these are in 
bijection with planted plane trees, and that under this bijection the height of a Dyck path corresponds to the height 
of the corresponding tree. The asymptotic behaviour of the average height of planted plane trees was determined 
by de Bruijn, Knuth and Rice [4], that is, they solved the average height problem for 1-watermelons with wall. 
Recently, Fulmek [9] extended their reasoning and determined the asymptotic behaviour of the average height of 
2-watermelons with wall. Katori, Izumi and Kobayashi [15] considered the diffusion scaling limit of 2-watermelons, 
and obtained the leading asymptotic term for all moments of the height distribution as well as a central limit 
theorem. The limiting process of p- watermelons has been investigated by Gillet [12]. He succeeded in proving the 
convergence of (properly scaled) watermelons to certain limiting processes, which he characterised by stochastic 
differential equations. 

In this paper we rigorously analyse the height of p- watermelons of length 2n with wall, and obtain asymptotics 
for all moments of the height distribution as n — > oo as well as a central limit theorem. In particular, we show that 
the s-th moment behaves like skI^tW 2 — S^K^j^n^" 1 )/ 2 + 0{n s / 2 ~ l + n p / 2 ~ p logn) asmoo for some explicit 

numbers «a , see Theorem [T] at the end of Section [31 The nature of our result explains the somewhat inconclusive 
predictions in [2]. To be more specific, Bonichon and Mosbah [2] predicted, based on numerical experiments, that 
Ki ~ \J (1.67p — 0.06). Although it does not seem unlikely that the constant k^ p \ as given in Theorem[TJ behaves 
like y/p as p — > oo, a rigorous proof is still lacking and work in progress. 

The proof of our result can be summarised as follows. As a first step, we represent the total number of 
watermelons and the number with height restriction in terms of certain determinants (see Lemma H]), the entries 
being sums of binomial coefficients. From these determinants we then obtain an exact expression for the s-th 
moment of the height distribution. After normalisation we may apply Stirling's formula and obtain an expression 
that can be asymptotically evaluated using Mellin transform techniques (see Lemma [7|). This kind of approach 
goes back to de Bruijn, Knuth and Rice [4]. Fulmek [9] adopted their approach for the asymptotic analysis of 2- 
watermelons. The new objects which arise here (and, in general, when extending this approach to the asymptotic 
analysis of p-watermelons) are certain multidimensional Dirichlet series (instead of Riemann's zeta function as 
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in [4] ) . An additional complication with which one has to cope is the increasing number of cancellations of leading 
asymptotic terms that one encounters in the calculations while the number p of branches becomes bigger. Thus, 
while a brute force approach will eventually produce a result for any fixed p (this is, in essence, what Fulmek [9] 
and Katori et al. [15] do for p = 2), the main difficulty that we have to overcome in order to arrive at an asymptotic 
result for arbitrary p is to trace the roots of these cancellations. We accomplish this with the help of Lemma [5] 
It allows us to exactly pin down which cancellations take place and to extract explicit formulas for the first two 
terms which survive the cancellations. The multidimensional Dirichlet series which arise in our analysis are the 
subject of the subsequent section. What we need is information on their poles. This information is obtained 
with the help of a relation that generalises the reciprocity law for Jacobi's theta functions (see Equation |3])), 
that is proved in Proposition [TJ We note that our definition of these Dirichlet series differs slightly from Fulmck's 
definition, which makes the analysis somewhat easier. These Dirichlet series that we encounter are related to 
so-called twisted multivariate zeta functions, studied, e.g., by de Crisenoy [5] and de Crisenoy and Essouabri [6]. 
However, their results cannot be used since they do not apply to our multidimensional Dirichlet series, which are 
explicitly excluded in these two papers. They can also be found within a class of multidimensional Dirichlet series 
studied by Cassou-Nogues [3]. In principle we could apply her results to our Dirichlet series and would obtain 
information on the poles of the analytic continuation of these series. But this would be cumbersome, and in our 
case it is more straightforward to obtain this information using the generalised reciprocity relation (see end of 
Section , which we are going to need in the proof of Theorem [T] anyway. 

We mention that small modifications immediately yield analogous results for p-watermelons with a horizontal 
wall positioned at some negative integer. Also, the analysis of the height distribution of watermelons without wall 
can be accomplished in a completely analogous fashion. 

The paper is organised as follows. The second section contains information on the analytic character of certain 
multidimensional Dirichlet series that is crucial for the proof of our main result. The third section contains the 
main result, see Theorem [T] at the end of that section. The techniques applied in that section are then used to 
obtain a central limit law, see Theorem [2] at the end of this paper. 

We close this section by fixing some notation. Vectors are denoted using bold face letters and are assumed to be p- 
dimensional row vectors. Further, we make use of the 1-norm and the 2-norm of vectors, viz. |w|i = wo + - ■ - + w p -i 
and |w|| = Wq + ■■■w p _ 1 . Finally, we define v w = Vq° . . . v p ^{ x . The relation v > w is to be understood 
component-wise. 



2. Some multidimensional Dirichlet series 
In this section we study the multidimensional Dirichlet series 



an O.p — 1 „ 

m n u ...m„_ 1 x - m a 



Za (z) = y .. 7 ^ , = E 



K + ... + m 2 y ^ | m |2»' 

where m = (mo, • ■ ■ , m p _i) ranges over IP \ {0}, for a = (ao, . . . , a p -x) € V ', a > 0. Our goal is to establish 
the analytic continuation of Z a (z) to a meromorphic function and the determination of its poles. Also, we need 
information on the growth of Z a (z) as \z\ — * oo in some vertical strip. 

It follows from the definition that Z ao __^ apl (z) = ^ a „ (0) ,...,o (T(p _ 1) (z) for every permutation a £ S p . If p = 1 
then 

Z a (z) = 2 [a even] £(2z — a), 
where [Statement] is Iverson's notation, that is 



[Statement] 



1 if 'Statement' is true, 
otherwise. 



If a p -i is odd, the definition shows that Z ao ^_ >ap 2tap l {z) — 0. Consequently, we may assume that the parameters 
ao, . . • , a p _i are even. 
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The analytic continuation of Zi & (z) is accomplished very much in the spirit of one of Riemann's methods for 
C,{z) (see, e.g., [25, Section 2.6]). In fact we have 

(1) {2m) T W Z 2a (z)T(z)= fV 1 I (n^W) -[a = 0])dfc, 



where $ a (t) = a (0, it) and where 



6 a (x,y) = —6(x,y) = £ {2mn) a e 2 ^ n+n *' 2 \ > 0, 

n— — oo 

is the a-th derivative with respect to x of y) = ^2 n e 2m ( xn + n2 v/ 2 ) ^ a variant of one of Jacobi's theta functions. 
Here, Equation ^ is obtained by substitution of Euler's integral for the gamma function, viz. T(z) = t z ~ 1 e~ t dt, 
and the series definition for on the left hand side of the equation above followed by interchanging summation 
and integration as well as a change of variables in the integral. 

We are now going to extract information on the poles of Zi & {z) from the integral (fTJ). This task is accomplished 
with the help of a generalised reciprocity relation (see Corollary [1]) , which is a consequence of the following two 
results, stated in Lemma [T] and Proposition [T] This relation generalises Jacobi's reciprocity law for 6(x,y), and is 
proved following along the lines of the proof of the reciprocity law in [19, Section 2.3]. 

Lemma 1. Let (f a (x,y)) a >o be a sequence of functions which are entire with respect to x for every fixed y with 
> 0. If {f a (x,y))a>o satisfies the conditions 

(i) f a (x + l,y) = f a (x,y) 

(ii) f a (x-y,y) = e^-y/V ± ($f k (x,y) 

then we have 



fc— n \ / 



cf\y) 



fc=0 

where 

pi 

c o k \v) = / fk(x,y)dx 
Jo 

is the constant term in the Fourier expansion of fk{x, y) as a function in x. 

Proof. Condition (0 implies the convergent Fourier expansion (f a (x, y) being understood as a function of x) 

fa(x,y)=Y, C n } (yy MXn+n2y/2) 
n 

for a > which shows that 

e -^- 3 //2) /a(x _ y)y) = ^ c W i( j /)e 2^n+nV 2 ). 

n 

Now, this last equation and Condition ([nj together imply the recursion 
which yields 



.k 



This proves the lemma. □ 
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Proposition 1. We have 

Proof. We prove the claim by applying Lemma [T] to the functions 

/.(*>») = E (--I'e^, a > 0. 

n ^ " / 

Condition (JTJ) of Lemma [T] is clearly satisfied by f a (x,y). For Condition we calculate 

f a{x _ y) = Wl _ — ] Q e — (-+»-«) 2 /y = e 2-(.-y/2) £ M £ f _£±^ ' 

n ^ » / fe=0 v / n v y 

It remains to determine the coefficients (y) of Lemma [TJ Short calculations show that 

4 a) (y) = jf f a (x,y)dx = 2 [a even] J ( -J e^^di. 
In particular we have for a = 

4%) = /"e- to, ^ = i 



g — i7r(x+n) /t/ 



Note that the evaluation of the integral above is true for y = it for some t > and analytic continuation then 
proves the correctness for general y with > 0. If a > then integration by parts yields the recursion 



and we obtain 

if a is odd 

¥ if a is even. 



4 a) (y) = 

Hence, by Lemma [I] we have 



(47riy)»/ 2 (a/2)! V i 



^ fa\(2k)\ j& (y\- k+l/ \ ( , 



On the other hand, expanding the binomial term shows that 



- kg - 2iri{xn-^n /y) 



k=0 

The last two representations for f a (x, y) prove the lemma. □ 
Putting a = in Equation ([3]), we obtain the reciprocity law for Jacobi's theta functions in the form 

J%e{ x ,y) = e-^'ve ( x -,-±). 

V * \V VJ 

Corollary 1. The functions r d a {y) = S a (0,iy), a > 0, satisfy the relation 

Proof. The corollary follows from Equation ([3]) upon setting x = and replacing y by i/y. □ 
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We can now prove the main result of this section. 

Lemma 2. The function Z^Jyz) can be analytically continued to a meromorphic function having a single pole of 
order 1 at z = | + |a|i with residue 



(5) 



TTP/2 

Res Z 2a (z) = —- — - T 

»=S+Wi r(§ + |a|i) 



n 

\i=0 



(2aQ! 
4 a >aA 



Furthermore, we have the representation 



[a = Ofr' , TT^IaU nUia 



Z2a(z) = - 



2 a i 



r(z + i) r(z) z 



Hi M)i«ur(z) A 



i*" 1 [ [ nW*)| -[a = 0])dt 



2-2|a|i 



(-4)i»iir(«) A 



*— M in*** 



/p-1 



_ (_^)|a| 1 ( TJ <&f. ) t^+^ ) (ft, 



where the two integrals above define entire functions with respect to z. For any non-negative integer k we have 

— 1 if a = and k = 0, 
otherwise. 



Z 2a (-fc) 



Proof. Consider again Equation JTJ), viz 

(2ttz) 2 I 



Z 2a (z)r(z) 



/°° ( fn ^ w j - 1 & = °i j *• 



We split the integral above into two parts, one over [0, 1] and one over [1, oo). The second integral is seen to define 
an entire function with respect to z. We consider the first integral 

jf ^ ( (ii ^ w j - t a = °i j d< = -^r^ + jf <z_1 fn w ) d< - 



By virtue of (HJ we obtain 



'p-i 



'p-i 



-p/2-|a|i 



0=0 



0=0 



= (-!)' 



<k<a \j=0 / 

+ (-^) |ai1 f n ^ ] ^ /2 - |a|1 ^ (-Y - 1 



, 0<k<a 



0=0 



Now, since for a^O the integrals 



1 dt 



define entire functions with respect to z we see that 



- /'''-'(n^G)) 



t -p/2-|a|i di 



defines an entire function with respect to z, too. 
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Combining all the parts and noting that 

/p-i/O m\ /•! (-_ N|a|! TTP-1 (2a,)! 

we obtain the representation for claimed in the lemma. The evaluations at the non-positive integers immedi- 
ately follow from this representation. □ 

We close this section with a result on the growth of ^2a(c + it) as \t\ — > oo. 

Lemma 3. For a e K fixed we have the estimate 

(6) Z 2 a(<r + ^) = 0(e £ l*l), |t|-K», 

for any e > 0. 

Proof. Mcllin transform asymptotics show that 

W)- ^a+1/2 ja+1/2 + ° I* )' * ~* ' 

^2a(t) = [a-0] + O(t- M ), t^oo, 
for any M > 0. Consequently, we have for a € W the asymptotics 

(fi (0) - ia - 0] - ^: n '^ +1/2> - [a - 0, + O (,») , ^ 0, 

/P-1 \ 



U^aM )-[a = 0}=O(t- M ), t^oo, 



\i=0 / 

for any M > 0. Now, by [8, Prop. 5] we see that the Mellin transform of (nf=ci ^2a(*)) — [ a = 0], viz- 

- 1 ^ a (^)r(z), 

satisfies 

/ 2 * a (a + ii) = 0(e-^ 2 - £ )l t l), |t|^oo 

for any e > and <r in any closed subinterval of (|a|i + p/2,oo), which can be extended to any closed subinterval 
of (—00,00) (see the proof of [8, Prop. 4] for details). The result now follows from the behaviour of the gamma 
function along vertical lines, viz. 

I> + it) ~ V2^\t\ a - 1/2 e-^ tl/2 , \t\ -» 00. 

□ 

3. The moments of the height distribution 

We denote by M^ h the number of p-watermelons with wall with length 2n and height strictly smaller than 
h. Further, we write for the total number of p-watermelons with length 2n. Note that = M^ h for 

h > n + 2p - 1 and M^ h = for h< 2p. 

Now, let Wn^ denote the set of p- watermelons of length 2n, and let P denote the uniform probability measures 
on these sets, and let H„ tP denote the random variable "height" on the probability space 

The goal of this section is to obtain an asymptotic expression for the s-th moment Ei?* p , where E denotes the 
expectation with respect to P, of this random variable as the length of the watermelons tends to infinity. Clearly, 
we have 

(7) Eff^ = _L_ £ (fc* - (fc - ( M &> - M$ h ) , s > 1. 

M 2n h>l 
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For determining the asymptotics of Eff* _ we proceed as follows. First, we find expressions in terms of deter- 

minants for the quantities h and M)f n . This is accomplished by an application of a theorem by Lindstrom- 
Gessel-Viennot, respectively of a theorem by Gessel and Zeilberger. Second, we obtain asymptotics for 

(8) M« fc and £ {h' - (h - !)•) (j|f&> - M^ h ) . 

h>l 

The asymptotics for EH^ p are then easily established. The main result is stated in Theorem Q] at the end of this 
section. 

(v) (v) 

We start with exact expressions for M^ l ' h and . 
Lemma 4. We have 

0) Mg= det (( 2U )-( . 

o<i,j<p \ \n + 1 - jj \n-l-i-j 

and 

(10) < h = Q Je^ ( ( n + ^ 1^ + l _ ] )-( n + m{h + 2 i r ;_ l _ l _ J ))), * > 0. 

Proof (Sketch). For h> 2p both equations follow from a theorem by Lindstrom-Gessel-Viennot (see [10, Corollary 
3] or [20, Lemma 1]), respectively from a theorem of Gessel and Zeilberger [11]. To be more specific, Equation ^ 
follows from the type C p case of the main theorem in [11], while Equation (fT0|) follows from the type C v case. 

The reader should observe that the entries of the determinant §§§ are the numbers of lattice paths from (0, 2i) 
to (2n, 2j) that do not cross the ir-axis. On the other hand, the entries of the determinant (fTUf are the numbers 
of lattice paths from (0, 2i) to (2n, 2j) that do not cross the x-axis and have height smaller than h. These sums 
are obtained by a repeated reflection principle (see, e.g., Mohanty [21, p. 6]). 

For < h < 2p the identity 



E 



2n \ / 2n 

n + rn(h + 1) + i — j J \n + m(h + 1) — 1 — i — j 



■E 



2n \ ( 2n 

n + m(h + l) + (h-i)- j) \n + m(h + 1) - 1 - (h - i) - j 



shows that the right hand side of (flO|) is equal to zero, since for h = 2i the i-th. row of the determinant is equal to 
zero, and for h = 2i+ 1 we see that the i-th and (i + l)-th row of the determinant only differ by sign and thus are 
linear dependent. □ 

We now turn towards the problem of determining asymptotics for the expressions ([5]). Asymptotics for the total 
number of watermelons are easily established since the determinant in ^ admits a simple closed form. The result 
is stated in the following lemma. 

Lemma 5. We have 




asm oo. 

Proof. The determinant © can be evaluated in closed form, e.g., by means of [16, Theorem 30], and is in fact 
given by 

M_ P t1 ( 2n t 2j ) f^Y (tt\o, . ^ fir 1 (2n + 2i)...(2n+l) 



m% = n 




o ( n+ n +1 ) V n J \fJ- ' ' ) \ f = L (n + 2j + 1) (n + 2j ) 2 . . . (n + l) 2 



This proves, upon determining asymptotics for the right-most product, the result as stated in the lemma. 

For a comprehensive discussion and references of this counting problem we refer to [18, Section 4]. □ 
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Asymptotics for the second part of |(HJ) are much harder to obtain. As a first step we note that 
(11) ^(h'-(h-W(Mg-M<*] h ) 



h>l 



= -Y(h'-(h-l)')Y det (( J n s )-( )) 

^ ^ Q o<i,j<v V \n + m t (h + l) + i- jj \n + m^h +l)-l-i-jJJ 

by an d (|10p. where the inner sum ranges over IP \ {0}. 

For determining asymptotics for (fTT|) we closely follow the proof of de Bruijn, Knuth and Rice [4] (in our case 
we have to overcome some additional difficulties). For sake of convenience, we give a short plan of the proof. As 
a first step we factor ( 2 ^) out of each row of the determinant on the right-hand side of (fl"U|) . We then replace the 
quotients of binomial coefficients by its (sufficiently accurate) asymptotic series expansion, which is determined 
with the help of Stirling's asymptotic series for the factorials (see Lemma E]). This shows that the asymptotic series 
expansion for (fTT|) can be expressed in terms of products of derivatives of Jacobi's theta functions we considered 
in the last section. With the help of the Mellin transforms and the results of the last section we are able to derive 
asymptotics for these functions (see Lemma [7]) • In Lemma [5] we finally obtain the desired asymptotics for (jlip . 

We start with the asymptotic series expansion for the quotients of binomial coefficients mentioned above. 

Lemma 6. For |m — z\ < n 5 ' s and N > 1 we have the asymptotic expansion 

I In \ 
/ir)\ \n+m — zl — m 2 /n 

{ ] ( 2 ;) 

/4JV+1 / x u / / \ 3Af+l u-1 21 / o \ /\/\ 2r+fc-l^ 



as n — ► oo. Here, the F r j are some constants the explicit form of which is of no importance in the sequel, and 
(— l) k k\<fik{w) is the fc-th Hermite polynomial, that is 

(13) Mz)= Y^(, m ){2z?™-\ fc>0. 

* — ' to! \k — ml 

m>0 v 7 

Proof. For sake of convenience, set x = (to — z)/n. With the help of Stirling's asymptotic series for the factorials 
we see that for x sufficiently small, \x\ < |, say, we have 

% -jrj- = [n + - j log — - nx log _ 

N , j l — 2k 

+ E ^Srr) ( 2 - f 1 + x ^ 2k - (1 - ^" 2fe ) + ( n "" 2iV ) 

for all fixed N > as n — > oo. Here, E>k denotes the fc-th Bernoulli number defined via z~2 k>0 B^t k Jk\ — tj (e* — 1). 
For the range \x\ < n -1 / 4 we further obtain by Taylor series expansion and some simplifications the expression 

2n \ /4N+3 yr \ i /4N+1 2 r\ 



, \n+m-z) I x r 1,1/ X T 



r 
r=l 



-1-2JV> 



(4W-1 / AT 

Further restricting ourselves to the range \x\ < n 3 / 8 we obtain, upon taking the exponential of both sides of 
the expression above and another Taylor series expansion, the asymptotic series expansion 

I 2n \ / 4N+3 /2AT-[3r/4J 

Vn+m — zJ 



-nx~ 



< |1+ Y I E l ]x 2r +0[7l 

i=-Lr/2j 



-l-2iV\ 
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for some constants F r j. 

Now, if N > 1 we obtain upon interchanging the two sums on the right-hand side above, replacing x with its 
defining expression (m — z)/n and simple rearrangements the expression 

2n \ / 3N+1 21 , N 2r 

e~ 



H^r = er{m ~ zfln ( 1 + E ^'EM^l + 

Vn/ \ 1=1 r=l 

Finally, expanding e~( m_z ) l n in the expression above in the form 



e -(m-z) /n _ e -m / r, 



•5> (£)(-£)'■ 



k>0 

and collecting powers of z, we obtain the result. Here, the </'fe(m/y / n) represent certain polynomials the explicit 
form of which is given in the lemma. □ 

We mentioned before, that the non-normalised s-th moment is a linear combination of certain functions 
related to products of the functions z?2a(i), a > 0, considered in the last section. In the next lemma we obtain 
asymptotics for these functions with the help of the Mellin transform and the results proved in the last section. 



2a 

m 



Lemma 7. For a 6 Z p , a > 0, and 4gN define the function 

(14) 5fe , a (n) = 5> + l) fc e ~ |m| = /n 

h>l m£(7i+l)Z p 
m#0 

For any fixed M > we have the asymptotics 

(15) 9kAn) = ( if ) ^ W + ^,a« (fc+1)/2 + (l - B k+1 t^j [a = 0] + O (n~ M ) 
as 7i — ► oo , where 

7— 1 + log-y^ ifp = fc + l, 
C(p — k) — 1 else, 

I j lim IZ 2a (z+|a|0r(z+|a| 1 )-(n 1 ^]f^) ifp=fc+l, 




and 



2 

' ' "T ^2a 

Here, 7 = 0.5772 ... is the Euler-Mascheroni constant 



rW + laMZaaf^ + lalO else. 



Proof. First, note that the function gk,a.(n) can be written in terms of derivatives of theta functions, viz. 

ffM( „) = (-^)-Ku S( A + i)*(M^) Wl (n^ (M!]_ [a = 0] 

h>i v / y=o v / 

Now, by the harmonic sum rule and Equation |T]), t ne Mellin transform of gk,a{n) is seen to be 

f°° 1 
/ g k< *( x dx = (C(2z-fe)-l)r(z + |a|i)Z 2 a(z + |a|i), > - max{p, fc + 1} . 

Consequently, the function g k ,a(n) can be represented with the help of the inverse Mellin transform by the contour 
integral 

9kA n ) = 7T~ / 5fc a (z)n z dz, c > - max {p, fc + 1} . 

Asymptotics are now being obtained by pushing the line of integration to the left and taking into account the 
residues of the poles of the integrand. 



9kA z ) 
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From the well-known analytic behaviour of the gamma and the zeta function (see, e.g., [26]) and the analytic 
behaviour of Zi & (z) as given by Lemma [5] we infer that the integrand g\. a (z)n z has potential poles at z = p/2, 
z = (k + l)/2 and z = — |a|i — m for m G N. For p =/= k + 1 all poles are of order one. Furthermore, the residues 
are given by 



Resgt & (zW- = (Cip-k) 1) I J] ||£ ] W /2 

\i=o 3 ) 

_§?!,/.*•«"■ = 5 r + Wl ) Zs - + |a| ') ••" + ' ),2 

Res 9^ a (»n z = - f ^— ^- - 1 ) [a = and m = 0] , 

\ fc+1 / 

where denotes the Z-th Bernoulli number defined via J2i>o Bit /V. — t/(e t — 1). 

In the case p = k + 1, the only difference is the pole at z = p/2, which is now a pole of order two. By Lemma|5] 
we know that 



Res Z 2a (z + |a| 1 )r(z + |a|0 = f JJ g^J I ,t" * 



and consequently, we have 



to, flSU(*)»- = fn ftrj (7 - l + logV^) (™) p/2 



Note that the limit above is equal to the constant term in the Laurent expansion of Z2 a (z + |a|i)r(z+ |a|i) around 
its pole z = p/2. 

For completing the proof we have to show the admissibility of the displacement of the contour of integration 
above. But this follows by well known estimates for the gamma and the zeta function along vertical lines in the 
complex plane together with Lemma [3] See [8] for details. □ 

This last lemma finally enables us to determine the asymptotics for the non normalised s-th moment (|lip . 
Lemma 8. We have the asymptotics 



h>l 

f 2n 



Jn~P 2 ^A i n'/ 2 -3^A a _in<"- 1 >/ 2 -|2P a (j[(2i + 1)M + O (i^' 2 - 1 + nP' 2 ^ 2 logra) 



as n — ► oo , where 



a>0 



with Wfc_i ja being defined in Lemma[7] 
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Proof. Substituting the determinant expressions and (fTU|) for and M^ h we see that 



^2(h s -(h-iy)(M^ 

h>l 

= Y(( h -iy-h°)Y det (( J 7 \ )-( ( 2 \ )). 

^ o<m<p \\n + m t (h + l)+i- jj \n + mi(h + l)-l-i-jJJ 

Instead of determining asymptotics for the right-hand side expression above directly we consider the more general 
quantity 

D n (x,y,z) = Y ((h-l) s ~h s ) Y det (Y 2 " )-( 2n )) . 

o<i,j<p \\n + mi + Xi— Vj / \n + m* — z — xt — w, / / 

Now, we factor ( ") out of each row of the determinant above, and restrict the sum above to those {p + l)-tuples 
(h, mo, . . . , m p _i) such that for i = 0, . . . ,p— 1 we have |(/i+l)mj| < n 1 / 2+e for some fixed e satisfying < s < 1/8. 
Since, by Stirling's formula, we have 



f 2 ™ ) 

= o 



whenever \a\ > n 1 ^ 2+E , we see that the sum of all terms failing to satisfy the condition above is O (n~ A/ ) for all 
M > and, therefore, is negligible. 

In the remaining sum we replace all quotients of binomial coefficients by their asymptotic series expansion 
as given by Lemma El Having done so, we extend the range of summation to Nx (V — {0}). This adds some 
additional terms, their sum being exponentially small and, therefore, again negligible. This technique of truncating 
the (exponentially small) tail of the exact sum, replacing the addends by their asymptotic expansion and finally 
adding a new (exponentially small) tail to the resulting sum has also been applied by de Bruijn, Knuth and Rice [4]. 

This procedure yields, upon noticing some cancellations due to summation over m which eliminates all odd 
powers of rrii for i = 0, . . . ,p — 1, for arbitrary N > the expression 



(16) Pat(x,v,z) 



h>l mS(/i+l)Z p 



m, 



where 

(17) T u , N {w,n) = cj> u {w) + Y ^EE^Iu-J^^ 1 " 



3JV+1 K-l 21 , _ 

2r 



such that 



(18) A l (x,y,z)=Q (p JV (x,y,z)+O(n- 2Ar - 1 G s , (n))), 

Here, the functions (j>k(w) are defined by (fT3|) . and 



. 2a 

m 



e -|m| 2 /n_ 



(19) G Sja (n) = £((fc-ir-fc s ) Y 

m^O 

Clearly, P/v(x, y, z) is a polynomial with respect to the variables xq, . . . , x p -i, yo, ■ • ■ , Vp—i, £• Furthermore, ex- 
panding the determinants and interchanging summations in (|16p reveals that P/v(x, y, z) is of the form 

(20) JMW)= 2 n (iwl|i|0/a (EgMJ^"" 1 ) '.^ 

i,j,l>0 \a>0 
|i|i + |j|i + |l|i even 
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for some polynomials <? a ,ij,i("- x ) in n 1 . Noting that [w 2a ]T2 U +i-N(w, n) — for all a by (fTT|) we obtain upon 
extracting the corresponding coefficients in (|16p the explicit representation 

(21) fe.v.w.^n- 1 ) = det ( ( Vl + Wj + U ) ((-l) Vi [k = 0] - l) K a *] T^.+^^n)) , 

0<ij<p\\ Vi,Wj,li J\ i J 

and short calculations show that 
[w 2a ] T 2u , N {w,n) 

= (-ir+^(u + a\ 8 ^ * ^ (-ij^+o-r / u + a - r \ 2o+2 „_ 2r _ fc 

(u + a)! V 2a 7 Z^ Z^ r > ! Z^ ( u + a _ r )| U a + 2u - 2r - k) 

v ' K ' 1=1 r=l fe=0 v ' v 7 

By expanding (/i — l) s — h s in powers of (h + 1) in (I19p and interchanging summations we see that 



(22) G., a (n) = ^ Q (2 s ~ fc - l)(-l) s - fc 9fc . a (n) 

where the functions gk,a{n) are defined in Lemma [7J 
Thus, we are led to consider sums of the form 



^ ga,v,w,l(n 1 )5fc,a(^)- 
a>0 

Now, we replace gk,&{n) by its asymptotic expansion (|15[) . viz. 

= fjl 1^4 j OfcW + - fe ,a- (fc+1)/2 + [a = 0] (l - 5 fe+ i^) + 

as n — > oo for all M > 0. The quantities f2fe(n) and Wfe ja have already been defined in Lemma [71 
The multi-linearity of the determinant in (|2ip then shows that 



vJ=0 



Ofe(n) det 



Wj + 10,' + / 



0<ij<p \ \ Vi,Wj,k 

The sum inside the determinant is further seen to be 

(2a)! 
4 Q a! 



(fo = 0](-ir h 2ai ] T^^^n) I 

a>0 ' / 



a>0 

3JV+1 2/ 2u-l 



f„-lV P V (-l) a+ - r (2r + fc - 2u)! v /k + r-u\( 2a 



2 2r+fc-2u (fc + r _ u )l Z-^ 7 ^ a J\2r + k-2u 

By the Chu-Vandermonde summation formula we obtain for the innermost sum above 
/ k + r — u\ ( 2a 



a>0 



2 



-Lfc/2J,^ + r-«+ Tfc/21 

i + \k/2] [k/2\ 



L(i + rfc/2l-Lfc/2j) r(«-r) 



r(i + rfc/2]) r(«-r-Lfc/2J) 

and, from the fact that [k/2\ > it — r on the right-hand side of the second to last equation above, we conclude 
that all terms having r < u vanish, since in these cases this last sum evaluates to zero. But this shows that 

(2a)' 



u/2 



n — > oo, 
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and we infer that 



a>0 



and further, noting that f2fc(n) = 0(n p / 2 logn) as n — + oo, 

(23) £ ?a,v,w,l(» -1 )flfc,«(») = E w fc,a9a,v,w,l(n- 1 ) n< fc+1 )/ 2 

a>0 \a>0 y 

+ go.v.w.iK 1 ) (l - B k+1 t^) + O („(i-Mi-lw|i-|i|0/a i ogn ) 



asm oo. 

Now, lets turn back to Equation (|16[) . Since the determinants involved in the definition of P/v(x, y, z) vanish 
whenever Xj = sGj or yi = yj for some i ^ j or i, = — z — xj or ?/, = — z — yj for some i and j, we conclude that 
Pat(x, y, z), which is a polynomial with respect to the variables x, y and z, is divisible by 

n & - x j)(yi - ) n ( z +^+^)( z +y i +2/j) 

^0<i<j<p J \0<i<j<p 

Consequently, all monomials of P/v(x, y, z) have total degree > 2p 2 . Furthermore, we see that 

(24) PAr(x,y,z) 

JJ [| (« + a* + aJj)(* + j n- pS C(n)(l + 0(n- 1 )) 

\0<i<j<p J \0<i<j<p J 

for some unknown function C(n) as n — ► 00. This function C(n) can be determined by comparing the coefficient 
of nf=o X T +1 vf +1 in © and (Ull). In this way we obtain 

(25) n~P 2 E «a,J,J I o(n- 1 )G. ) a(n) = 4%-P 2 C(n) (l + 0(n- v )) , 

a>0 

where J = (1, 3, . . . , 2p — 1). Since 



n — > 00, 



(26) G Sja (n) = -sg s _i )& (n) + 3( ]5 s _2, a (^) + 0(g s - 3 , a (ri)), n — > 00, 



by |22|). we see by (|23J) that 

E 9a,j,j,o(n- 1 )G s , a (n) = -£ ga.J.J.o^ 1 ) (sg a -i. a (n) - 3 ff) fl,_ 2 ,a(n) ) + O [n 8 ^ 1 + nV/ 2 -P 2 log7 

a>0 a>0 ^ ^ ' 

as n — > oo. Noting that 

2?(-4)Hi / (2i + 2j + 2)! , 

9a,J,J,o(n- ) = „ det + 0(n ), n -> oo, 

iT-o ( 2 * + 1)' ' - IJ<P V ( J + J + 1 - a 0! (2a l )!/ 



we further see by |23J that 



£ ?«,J,j,o(n- 1 )G., a (n) = ^ j ( sA ^ s/2 ~ 3^) A,-i^- 1)/2 + A + O (V 2 - 1 + n p / 2 - p2 logn)) 



as n — > 00, where 



A, = - E(- 4 ) |a|1 n< det ( r ^ (2 lt 2j+ vl' J * > 0, 

^ o<i,j< P + j + 1 - Oj)! {2ai)\J 



a>0 
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3 f(2i + 2j + 2)\ 

Aq = — — det 



and 

2 0<i,i< P V + 1)! 

Here, the constant Ao is of interest only in the case s = 1 (it can be absorbed into the O-term otherwise), and 
comes from the asymptotic expansion of gi : o (n). 

Now, with the help of Equation (|25p we can determine asymptotics for the function C(n), which gives us 
asymptotics for P/v(x, y, z) by Equation (f2~l|) . and finally also asymptotics for Z?at(x, y, z) by Equation (fT%|) . 

The proof is now completed upon specialising to a;, = yi — i for i = 0, . . . ,p — 1 and z = 1 in the asymptotics 
for Dn(x, y, z). For sake of convenience we finally note the identities 

n ( j -j) 2 ) ( n (i+^+j) 2 ) =n(2 l+ i)! 2 , 



det '( 2i + ^ + 2 ) ! 



.... . . ?^)=2^n(2i + l)!. 

The second identity can be proved by means of standard determinant evaluation techniques (see [16] for details). □ 
Finally, we can state and prove the main result of this paper. 

Theorem 1. Set M p = 2 p2 l\iZo( 2i + x ) ! and T pW = det ($2i+2j+2(t)). For s e N, the s-th moment of the 

o<j,j<p 

height distribution of p- watermelons with wall satisfies 

(27) EHl p = SK Wn s / 2 - 3 Q K^-D/a _ | + ( n ./a-i + n p/2- P 2 bg ^ 

as ?^ — ► oo , where 

Proof. Replacing and the sum in Equation with their asymptotic expansions as given by Lemma [5] and 

Lemma |5] we see that 



Q^-V*" 1 ^ - f + O {n s / 2 - 1 + TiP/ 2 ^ 2 logrc) 



where, for k > 0, 



K (p)-_^Vf-4)l^ det f (2i + 2j + 2)l \ 



The quantity Uk-i^a has already been defined in Lemma[7J 

In order to prove the integral representation for when s ^p, where we have u>k-i,a = |T (| + |a| i) i?2a (f + |a|i) , 
we consider the more general expression 

(*) =-^rE n < det f (2 : + ^ + 2)!( 7f;, ) r (« + |a|0 Z 2a ( Z + \a\i) 
2M B z — ' o<ij< p V (z + 7 + 1 - a. ! (2a,; ! / 



tp^ Q o<i,3<P \ + j + a l )\(2a i y. 

(2z + 2j + 2)!(i/7r) a 



2^pJo \ ^(t+j + l-o)!(2o)! 



l?2o(*) - M p ] dt 



for 5ftz sufficiently large. Here, the second line is a direct consequence of Equation ([T]). The reciprocity relation ^ 
followed by the change of variables t^t^ 1 then shows that 
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Table 1. This table gives numerical approximations for sn f for small values of s and p. The 

(v) 

quantity sk; s is the coefficient of the dominant part of the asymptotics for the s-th moment of 
the height of p- watermelons (see Theorem [T]). The calculations have been carried out using the 
integral representation for as given in Theorem [T] The results shown here conform well with 
with numerical results obtained by Fulmek [9] and Katori et al. [15]. 





s = 1 


s = 2 


s = 3 


s = 4 


P= 1 




3.289... 


6.391... 


12.987... 


p = 2 


2.577... 


6.790 . . . 


18.282... 


50.306... 


p = 3 


3.207... 


10.429 . . . 


34.371 . . . 


114.817... 


p = 4 


3.742 . . . 


14.141 . . . 


53.939... 


207.712... 


p = 5 


4.215... 


17.898... 


76.536 . . . 


329.655... 



Asymptotics for $2a(i) for t — > and t — > oo as given in the proof of Lemma [3] then show that this last integral 
is convergent for 5Rz > 0. The representation for s ^ p is now proved upon observing that, by definition, we have 

Now, consider the case s — p. Here, we have 



w p -i,a = lim \Z 2a .(z+ |a|i)r(z + |a|i) 

Z \ 

which shows that, as in the other case, 

In this last calculation, we have, after interchanging the order of the limit and the sum, applied the results obtained 
in the case s ^ p. This proves the theorem. □ 

Some numerical approximations for the coefficient of the dominant term of the asymptotics proved in Theorem[T] 
are shown in Table [3] But our last theorem does not only give the dominant term of the asymptotics of the s-th 
moment of the height distribution of p-watermelons but also the second order term. So, for example, we obtain 
the more precise asymptotics 

■ I I, A 

Remark 1. It can be shown that Theorem Q] is even valid for s G C, 5i(s) > 0. The proof of this more general 
result is the same as for our theorem except for two small changes which we are going to address now. 

In the proof of Lemma[8]we defined the functions G Sta (n) (see Equation (|19p). For s G N the asymptotics ([26]) 
for G StSL (n) were easily found by the expansion p2|) . This is not possible for s G C \ N. In order to prove the 
asymptotics (I2TJ1) in that case we note that (see Equation |H])) 

The term for h = 1 in (flT)]) is seen to be negligible due to summation over a > (see the discussion of the function 
f2fc(n) following Equation (f2"2"|) in the proof of Lemma [H]). For h > 2, we can use the binomial series expansion in 
the expression above and finally obtain the asymptotics (|26p . 

The second change concerns Lemma [71 which has to be generalised to k G C. But this makes no difficulties. 



** /a (ng m) \ 
(--I) J 



= \/nn — — + O yn log nj , n — > oo, 

= 2.577... ~ 7j + (™~ 1/2 ) , n-^oo, 
= 6.790... n- 3.866... Vn + 0(1), n->oo. 



THE HEIGHT OF WATERMELONS WITH WALL 



17 



4. A CENTRAL LIMIT LAW 

We are going to derive asymptotics for the cumulative distribution function of the random variable "height" on 
the set of p- watermelons with length 2n with wall, i.e., 

F n (h) IP { //,..,. - h] ' 



for the range h + 2 = t^/n, where t £ (0, oo) 



Theorem 2. For t € (0, oo) fixed, the random variable H n ^ p on the set of p-watermelons of length 2n with wall 
satisfies 

{ ' \ y/n - ) (~2)P 2 nr o 1 (2i + l)!o<M< P V 2l+2l+2 \t 2 JJ \nt 

as n — > oo, where the constant implied by the O-term is independent of t. 

Proof. The result can be proved in pretty much the same way as Theorem [T] Therefore, we only give a rather 
brief account of the proof, and refer to Lemma EI for the details. 

Instead of the exact expression (JTUJ) for M^ h we consider the more general quantity 

p I I 2n \ I 2n \ 

2fl \ j ^ | y "v \n-\-m-\-Xi— yj) \n+m—z—Xi—yj) 



(29) F n (h; x, y,z) — ( det V 

\ n I ft<7 t^rt \ ' 

v me(/i+2)Z 

Again, we find the polynomial 



n J o<iJ<P \ ^ ( 2n ) ( 2n ) 

2N 



Jv x,y,z )= det > J - > T 2u -n \-i=,n] e 

\u=0 mS(/t+2)Z 1 v 



such that 



F n (h; X ,y,z) = ( 2 "Y ( Q N ( X ,y,z)+0 ( n^ 1 £ • ' " 



e 

m£{h+2)Z 

where N can be chosen arbitrarily large and T2 U -n being defined by (|17p . The polynomial Qat(x, y, z) is seen to 
be divisible by 

n & -^fa -Vj) i n (s+^+^'X^+s/t+fj) 

v o<i<j<p y yo<i<j<p 

since the determinant in the definition of Qn(x, y, z) vanishes whenever Xi — Xj or g/j = y 3 for some i ^ j or 
Xi = — z — Xj or 2/j = —z — yj for some i and j. Hence, 

Qjv(x,y,z)= JJ (jBi-Sj-Jd/i-yj) J II {z + x l +x J ){z + y l + yMC{t)(l + 0{n- 1 )) 

\0<i<j<p J \0<i<j<p J 

as 7i — ► oo for some unknown constant C(t). Now, we are going to determine asymptotics for C(t) as n — > oo. 
This task can be accomplished by comparing the coefficients of the monomial Ilo<i<p X T + ^ Vi l+1 m ^ ne expression 
above and the defining expression for <5at(x, y, z). We obtain 

4PC(t) (1 + 0(n-i )) = (-2)^ 2 dot ( ( 2i + ^ + 2 ) £ T 2i+2j - +2;JV (mt, n) e~^A . 

Recalling the definition of the functions T2 a ;iv(w, n) (see Equation (|17p), our attention is drawn to sums of the 
form 

oo 

(2a)! 4>2a{mt)e-^ 2 , a e N, 

71,=— OO 

where the polynomials 02a are defined by (jTT 
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Now, rewriting the reciprocity relation (U|) as 



a(-) =y a+1/2 7T a J2 (2a)!0 2a (n0Fy)e-" 2 ^ 



we see that 

(2a)! £ ^(mt)e-(-) 2 = ^ 2a (J) . 

n— — 00 

From the asymptotics for $2a(£) as given in the proof of Lemma [3] we deduce that 

= const + 0(t-*), t-^oo, 

^ 2o g) = [a = o]f + o («-"), t-Q, 

for all M > 0. Consequently, we obtain for h + 2 = ty/n, where t is fixed, 

E W ("*• «) e_(mt)2 - ^TT^a (J) +o(±), n ^ co. 

Note that the constant implied by the O-term can be chosen independent of t. 

The theorem is now proved upon substituting these asymptotics for our sums appearing in the expression for 
C(t) above, taking out some factors, specialising to Xj = ijj = 1, and dividing by ■ D 

Remark 2. The asymptotic cumulative distribution function of the random variable "height" as given in The- 
orem [5] has been re-derived by two groups since the first version of this manuscript was distributed. Since their 
expressions differ from the one given by Equation (|28p we want to give some comments on the equivalence of these 
three (more or less) different expressions. 

The expression found by Katori et al. [14] can easily be obtained by an application of the reciprocity relation ((4]) 
to Equation (|28p . and therefore, is not essentially different from the one given here. 

Schehr et al. [24] expressed the cumulative distribution function of the height as a multiple sum, which can 
also be easily derived from Theorem [2] Since i?2a(^) = X)^L-oo( — 47r 2 n 2 )°e _n %t we see that the determinant in 
Equation is equivalent to 



det U 2i+2j+2 (J) ) = V e -^oW/*) 2 det ( { -4 n yy+^) 



no,...,ra„_i6 




( _ 7r2) P 2 2 2P 2 +p £ yi n/ +2 ] ( n k 2 - »?) ) e E ^ i(n ^ /t)2 

n ,...,n p _i£M \j=0 



l<no<-<rip_i \0<i<j'<p / \a£& p j=Q 



{-■k 2 Y2 2p +p 

l<no<"'<n p -i \0<i<j<P 




where & p denotes the set of permutations on the set {0, 1, ... ,p — 1}. Here, the second and fourth lines are simple 
consequences of the Vandcrmonde determinant formula. For the third line we note that the the factor sgn (a) 
cancels the sign change of the Vandermonde product due to the rearrangement from the second to the third line. 
Substituting this expression for the determinant in Equation (|28p we arrive at 

2 

OP 2 +P 7T 2p 2 +p/2 f -2p 2 -p ,' \ , „ 



n „-. (2j + 1) , E | n <»?-?>, in. 



l<no<-<™ P -l \0<i<j<p 

which is Schchr's expression for the cumulative distribution function 
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For the special case p = 1 we obtain the well known central limit law first proved by Renyi and Szekeres [23] , 
viz. 

(30) P {^T - 4 £ (l-2(m0 2 )e-( mt ) 2 , n^oo. 

It should be mentioned that this limiting distribution is the distribution function of v^maxo^^^i e(x), where e(x) 
denotes the standard Brownian excursion of duration 1 . For details and references we refer to the survey paper by 
Biane, Pitman and Yor [1], in which the authors consider probability laws related to Brownian motion, Ricmann's 
zeta function and Jacobi's theta functions. 
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